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Abstract 



Oh! 

^ ■ Proton decay via d=5 operators is excluded by now not only in the framework 

of SUSY SU(5) but also its extensions like SUSY SO(IO) are on the verge of being 
inconsistent with d=5 decays. It is reasonable therefore to suppress, e.g. by a sym- 
metry, the d=5 operators and to consider gauge boson induced d=6 decays. This is 
suggested in several recent papers in the framework of models with a lighter Mx ■ We 
C^ ' discuss here explicitly the fermionic sector of such a renormalizable SUSY SO(IO) with 

realistic mass matrices. We find that the recently "observed" large leptonic mixing 
leads to an enhancement of the nuclcon decay channels involving /i's and in particular 
the fJ-'^TT°, /i+TT^ modes. 



*e-mail:achiman@theorie. physik.uni-wuppertal.de 
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Nothing is known directly at present about the supersymmetric (SUSY) partners, except 
for the experimental lower bounds on their masses. The main indication for (low en- 
ergy) SUSY is the unification of the coupling constant of the minimal SUSY Standard 
Model (MSSM) at a high scale « 10^^ GeV. However, this unification speaks even more 
for SUSY-Grand Unified Theories (GUTs) |j|], which predict also the observed partial 
Yukawa unification (i.e. rrir = mi ) at those energies. Yet proton decay, the other main 
prediction of GUTs, was not observed till now [^. Is this consistent with the expectation 
from SUSY GUTs? 

It is generally assumed that the leading proton decay modes in SUSY-GUTs are due to 
d=5 operators, as the d=4 contributions must be suppressed via a symmetry (like R- 
parity) to avoid a much too fast decay. Those d=5 decay modes are induced using the, 
not yet observed, SUSY particles and involve therefore many unknown parameters [^-[0]. 
As a consequence the predicted rates are highly model dependent. This large freedom does 
not save SUSY-SU(5) from being practically ruled out if the needed threshold corrections 
are used to limit the mass of the Higgsinos {Ht) which mediate proton decay Q ||^]. Even 
extensions of SUSY SU(5) like SUSY SO(IO) are on the verge of being excluded by similar 
arguments |^]|^. To save it one must push several parameters to their allowed limit and in 
particular M^ is required in some papers to be much larger than Mqut and even larger 
than Mpianck = 2.4 x 10^^ GeV. This leads to very large couplings in the superpotential 
or the corresponding non-renormalizable contributions. 

In this letter we will assume, as was already put forward by several people P |0], to 
suppress not only the d=4 but also the d=5 contribution by a symmetry and consider the 
gauge mediated d=6 operators. 

• Those involve only known coupling constants and masses and therefore their pre- 
dicted rates are far more reliable than the d=5 ones. 

• They are the real test of the GUTs because d=5 induced proton decay is allowed in 
non-GUT theories as well. 

• Although d=6 contributions are suppressed by 1/M^ they could yet be observable 
in the near future if the relevant gauge bosons mass Mx is somewhat lower than 
10^6 GeV. 

This gauge-boson induced proton decay does not require non-renormalizable contributions 
or particles with masses above the GUT scale. We will apply therefore the "renormalized 
see-saw mechanism" using *I*y26 ^^ Sive large masses to the right-handed (RH) Majorana 
neutrinos. This has many advantages, in particular R-parity invariance, that is needed 
to avoid the catastrophic d=4 contributions, is automatically obeyed in this case. This 
invariance leads also to a stable neutralino as a natural candidate for dark matter. On top 
of that, a way to lower the GUT scale in terms of a fully renormalizable gauge theory was 
suggested recently by Aulakh, Bajc, Melfo, Rasin and Senjanovic |10|. They introduce 
an intermediate Pati-Salam gauge group [^ at a scale Mj. Such an intermediate scale 
is needed anyhow to explain the fact that the masses of the heavy Majorana neutrinos 
required to give the light neutrinos masses consistent with the observed oscillations fl^ . 



^See however Refs. 



are considerably smaller than the unification scaleFI. In this SUSY-SO(IO) model, due 
to the absence of trilinear terms in the superpotential, some particles acquire masses of 
0{M'j /Mgut) via mixing. Those particles affect the renormalization group equations in 
a way that lowers the unification mass. The model can solve the doublet-triplet problem 



a la Dimopoulos-Wilczek [15| and can suppress naturally the d=5 operators i.e. Higgsino 
mediated proton decay. 

Mx and/or the unification scale have smaller values also in other recent models Q and in 
particular in those using large extra dimensions. As a specially interesting example let us 



mention here the paper of Hall and Nomura |16| based on the model of Kawamura [17|. 
The idea is to use a five dimensional SU(5) GUT compactified to four dimensions on the 
orbifold S"^ /{Z2 x Z2). This yields the MSSM with doublet-triplet splitting and a vanish- 
ing proton decay from d=5 and d=4 operators by a U{1)r symmetry. The model gives 
a compactification scale {Mc = Mx) somewhat lower than the four dimensional unifica- 
tion scale. This model is not a 4d GUT, but it involves all the properties we use for the 
fermionic mass matrices and proton decay. 

When d=6 contributions are discussed in the literature one refers always to the proton 
decay into e^i:^ , that is the dominant decay mode only when the mixing is neglected. 
This is in contrast with recent papers about d=5 proton decay H H, where the effects of 
realistic mass matrices are explicitly considered. 

We present here a model for the explicit realization of the fermionic sector of a SUSY- 
SO(IO) GUT with realistic fermionic mass matrices. We will in particular show that the 



observed |14] large leptonic mixing leads to the enhancement of the branching ratios of 
the nucleon decay into muons with respect to those with e^ in the final state. 
The observed large leptonic mixing will be coupled in our model with large quark mixing. 
Note, that the CKM matrix gives the difference between the mixing angles of the up and 
down quarks. Therefore, only the difference between the lefthanded (LH) mixing angles 
of the quarks must be small while the righthanded (RH) rotations are unobservable in the 
framework of the SM or the MSSM and can be large. The nucleon decay in GUTs is one 
of the few observables in which all mixing angles are involved^. 

The aim of our effective model is not to get as many as possible "predictions" for the 
known observables of the SM. It is to calculate all the "non-observable" mixing angles in 
terms of the observables of the SM in order to predict the proton decay. 



The model uses a scheme developed in a series of papers [^ [£1| . All mass matrices have 
a non-symmetric Fritzsch texture ||2^ [23| 




M = I B C \ . (1) 



This texture and the contributions of the Higgs representations are fixed by a global Uf{1 
(or Zn). We use only renormalizable contribution a la Harvey, Ramond and Reiss [p4| 



^This scale is needed also for the invisible-axion ||l3] and the baryon asymmetry induced via the lepto- 



genesis TL 
^Partii 
leptoquark interactions, baryon asymmetry due to leptogenesis e.t.c 



^Partial information on the "non- observable" angles can be obtained by looking for RH currents, scalar 



This has some interesting advantages in our opinion, as was mentioned before, compared to 
the method of Froggatt and Nielsen [^] that uses broken Uf{^)^s with non-renormahzable 
contributions Q. 

The SO(IO) symmetry as weh as Uf{^) give relations between the different entries of the 
mass matrices. We need one large VEV in <&y26 ^° &^^ ^^^ ^^ Majorana neutrinos masses 



the corresponding mass matrix must have the texture |19| 




M,^ = a . (2) 



Because our main interest lies in the nucleon decay, CP violation is neglected in this letter 
and the parameters are taken to be real, for simplicity |^. 



We use, as is discussed in [^, one heavy ^y2E *o &^^ the RH neutrinos a mass. To gener- 
ate the light mass matrices one $j26 and pairs of $10 and 'I'j^o ^^^ needed f\ Those Higgs 
representations involve several SUl{2) doublets. Note, however that as is usually done in 
SUSY GUT models |l8| pl broken into the Minimal SUSY Standard Model(MSSM) at 
the TeV scale, only two combinations of those doublets remain "light" . All other com- 
binations acquire heavy masses and therefore do not affect the gauge unification |^. The 
explicit combinations are not important because the effective VEVs are anyhow free pa- 
rameters in our model. For the ratio of the light VEVs we use in this letter the commonly 
used value tan/3 = 5 ^ The gauge coupling becomes non-perturbative soon above the 
unification point but we do not use physics above the GUT-scale. 

The three fermionic families in 16i i=l,2,3 and the Higgs representations ^k transform 
under the global U{1)f as follows: 

16j — >■ ex.p{iaj6)16j (3) 

$fc -^ exp(i/3fc0)$fe. (4) 

Invariance of the Yukawa coupling terms 16i$fcl6j under U{1)f requires Oi + Oj = j3k- 
Hence, the most general structure of the Yukawa matrices is: 



Y« = XI ; Yi^,), = Ui ; Y« 






Y(^„' =0 X2 ; Yi^,', = ya ; Yg„ = Z2 • (5) 



These Yukawa matrices give explicit expressions for the M^, Mu, Me and MJ^ mass 
matrices, in terms of a set of 14 parameters (combinations of the Yukawa couplings and 

''The method of Froggatt and Nielsen is useful to explain the hierarchy in the quark mass matrices but 
not for the neutrinos. This is because the matrix elements are fixed in this method only up to unknown 
0(1) factors and the see-saw matrix is a product of three matrices. The neutrino matrix elements are 
given in this case only up to corrections of [0(1)]^ which can be large. 

^Complex parameters will be used in a forthcoming detailed paper pol] 

®The model involves pairs like 126 + 126 e.t.c. to avoid breaking of the SUSY at high energies, but 
only part of those are relevant for the fermionic mass matrices. 

^Some light VEVs are then induced ones. 

*In Ref. pel we shall discuss a set of diff'erent values of tan /3. 



the corresponding VEVs). These matrices are diagonalized and fitted to the the known 
masses and mixing of the quarks and leptons as wih be described in the fohowing. 

We start by taking diagonal mass matrices for the charged fermions at Mz- Than the 
fuh two-loop renormalization group equations (RGEs) of the Minimal Supersymmetric 
Standard Model (MSSM), with tan/3 = 5, are used to get the following masses at Mqut = 
2 X 10^^: 



m„(MGUT) 


"id(MGUT) 


m^(MGUT) 


1.04 MeV 


1.33 MeV 


26.5 MeV 


mc(MGUT) 


mb(MGUT) 


mi(MGUT) 


302 MeV 


1 GeV 


129 GeV 


me(AfGUT) 


m^(MGUT) 


m^(MGUT) 


0.32502032 MeV 


68.59813 MeV 


1171.4 MeV 



In our model only 12 independent parameters are needed to specify charged fermion ma- 
trices at the GUT scale, exactly the number of the underlying masses and CKM mixing 
angles (the CKM matrix changes only slightly from Mz to Mqut)- By fitting those pa- 
rameters, all the LH and RH mixing angles of the quarks and leptons are fixed[|. This 
procedure involves a set of non- linear equations: 



utM„U, = 


= M(f), 


dIm.d, = m'P 


e1m,e^ = 


= M(^) , 


UtD, = VcKM . 



(6) 



We found five solutions to those equations all of which have several large mixing angles [^0[ . 
These fix the neutrino mass matrices M^**" and Mj^j^ also up to two parameters (and the 
overall scale M/j). Hence, the see-saw matrix 



M' 



light 



-m'^^'^^ {M.'^^'^'^Y^{mf''^y 



(7) 



as well as the leptonic LH mixing \Jmns = E|_Nj^ are also known. We varied than 
the two free parameters and looked for solutions that reproduced the neutrino data with 
LMA-MSW or SMA-MSW for the solar neutrinos p|]. The details will be given in a 



forthcoming paper |20]. The best fit is obtained for the a LMA-MSW solution with the 
following properties: 



QU 

od 

od 
^R12 



^L12 



-1.1 X 10"^ 



1. The Quark LH and RH mixing angles at the GUT scale: 

ei,, = -0.077, ei,, = -1.48, ei,, = -4 x lo-^. 

-0.045, e-^,, = -2.2xW-\ 0«,3 = 
0.15, 0'[,^ = -l.U, ^^,3 = 1x10-5. 
-0.33, ^^23 = -3xl0-^ ^^13 = 6 X 

2. The mixing angles of the Charged Leptons: 

-1.17, eL. = 1.44, eU. = 0.0002. 

0.002, 



10 



-2 



'R12 



"L23 

9^ - 
^R23 ' 



'L13 

-0.003, 6l|j^3 = 0.002. 



If only one parameter is taken to be complex, one can use its phase to account for the observed CP 
violation and again all mixing angles will be fixed pffl. 



3. The Neutrino masses : 



Mr = 5.2 X 10^3 GeV 



m^^ = 1.88 X 10-3 eV, m^^ = 5.89 x lO^^ eV, m^^ = 5.85 x lO^^ ^y. 

4. The LH Leptonic (Neutrino) mixing angles: 
6*5^2 = 0.55, 61^3 = 0.74, 9^^ = -0.0053. 

Using these results one can calculate the proton and neutron decay branching ratios. We 



use the method of Gavela et al. ||2^ as it was extended in a series of papers pq| |21] to 
models with large fermionic mixing (especially RH ones). In this work it is once more 
generahzed into a SUSY GUT. 
The effective baryon number violating Lagrangian of SO(IO) is 

+ ^3 (£<.,,nf-7"0 ifitl.dl) + A^ (e„,,nf-7"0 (AJTmO 

+ A{e^,,u'i''r<){et7,sl) + A,{e^,,uf''r<){etj,s-^) 

+ A, {e^,,uf^r<) ifltl.si) + ^8 {e^p,u1-rO (A^7mSS) 

+ A9(e.,,nf-7"df)(i>f«7M^S) + A,o{e^p,u'l'7''di){i?^^j,d%) 

+ A,,{e^,,u^-Ydl){i>f^7,sl) + A,^{e^,,u1''rdt){iy^^7,s%) 

+ A,^{e^,,ut-rsl){i^^n7,d-^) + A,^{s^,,ut'7'si){i^^^j,dl) 

+ A,^{e^,,ut'rdl){i^^^^,dl) + A,^{e^,,ut'rdi){i^?^7,sl) 

+ A,, {e^,,ut-rsl) {i^?nl.d%) 
+ ( terms with two s quarks ) 
+ ( terms with c ,6 and t quarks ) 
+ ( terms with f +^ and pf;^^^ ) 

+ h.c. (8) 

(9) 

where Ai are functions of the mixing angles given in Ref . ^l| . 

The partial decay rate for a given process nucleon — > meson + antilepton is expressed as 

follows 



I r 

where M^ and Aij. are the hadronic transition matrix elements for the relevant decay 
process. I and r denote the chirality of the corresponding antilepton |28|. Ai and A^. are 



the relevant coefficients of the effective Lagrangian (N). A, A^ and Aj^ are factors which 



result from the renormalization of the four fermion operators, see Ref. [20| 



Using all this we obtain the nucleon decay branching ratios for our best solution. These 
are presented in the following table compared with the d=6 nucleon decays without mixing. 



proton 


% 


% 


neutron 


% 


% 


decay channel 


no mixing 


LA-MSW 


decay channel 


no mixing 


LA-MSW 


p -^ e+vr'^ 


33.6 


17.5 


n -^ e+7r~ 


62.86 


32.5 


p -^ /i+7r° 


- 


16.1 


n — > /i+TT^ 


- 


30.0 


p -^ e+K^ 


- 


4.6 


n -^ e+p^ 


9.7 


5.0 


p -^ n+K^ 


5.8 


2.7 


n -^ P^P^ 


- 


4.6 


p -^ e+r/ 


1.2 


0.6 


n -^ z/fTfO 


15.1 


9.2 


p -^ II+T] 


- 


0.6 


n -^ ufK'^ 


- 


2.6 


p -^ e+/9° 


5.1 


2.7 


n — > z/fr/ 


0.6 


0.3 


p ^ ^v° 


- 


2.5 


n ^ <7r0 


- 


5.1 


p -^ e+w 


16.9 


8.8 


n -^ ly^K'^ 


1.7 


0.0 


p -^ fl^UJ 


- 


8.1 


n -^ i^^rj 


- 


0.2 


p -^ i/f 7r+ 


32.3 


19.7 


n ^ z^fpO 


2.3 


1.4 


P -^ <7r+ 


- 


10.9 


n -^ u^LO 


7.7 


4.7 


P -^ <^+ 


0.1 


0.2 


n -^ u^p^ 


- 


0.8 


P -^ i^f P'^ 


4.9 


3.0 


n -^ v^LO 


- 


2.6 


P -^ y^K*+ 


- 


0.1 


n — > u'^TT^ 


- 


0.1 


P -^ v'ilP'^ 


- 


1.7 


n -^ u^K^ 


- 


0.7 



One sees clearly that the nucleon decay rates into muons are strongly enhanced. Also 

decays into e'^K^ and vK^ are not negligible. 

The other solutions in the LMA-MSW case give similar results. Since we use the most 

general Yukawa matrices (in terms of the asymmetric Pritzsch texture), we expect our 

results to be generic in the large mixing scenario. 

The absolute decay rates are much less reliable than the branching ratios given above. 

They depend not only on the "unknown" value of Mx, but also on the uncertain hadronic 

matrix elements p^], the short distance enhancement factors e.t.c. Using estimates for 

the absolute decay rates |^] ||7|, in terms of the recent Lattice calculations of the JLQCD 

collaboration pO] , we obtain for our branching ratios the following result 



1 



1 



i« in34 ( Mx ^V0.015GeVS2 
16 X 10 X I — TT^- I 



r(p ^ /i+7r°) " r(p ^ e+7r°) — " — "KlO^GQ^yJ 
One can apply now the Super-Kamiokande bound 0] 



an 



yrs. 



— — - > 4.4 X lO-^'^ yrs (90% CL) 

r(p -^ e+TT") ^ ^ ^ 

to get a lower bound on effective Mx's in the theory 



Mx > 0.7 X 10^^ GeV . 

To observe gauge mediated proton decay in the near future one needs surely lighter X- 
bosons than 10^^ Gev. This is a natural feature of several recent models as indicated 
before. Super-Kamiokande will continue to look for the nucleon decay and expects to 



reach in 8-years the limit of 1.4 x 10'^^ years |31]. Some of the new nucleon decay ex- 
periments and especially ICARUS l35] and the half megaton UNO detector [p3| are well 
suitable to look for decays into /i's and i^'s. Therefore, if the contributions of the d=5 
operators to the nucleon decay are really suppressed the search for the fi~^TT", /.i^vr^ modes 



should not be neglected. 

The enhancement of the muon branching ratios is a unique feature of our model because 
the decay mode p — > e^ir" is not negligible also in the d = 5 induced decays [Q]. In view 
of the fact that this enhancement is the effect of the large observed leptonic mixing on the 
d = 6 nucleon decay, we suggest that the observation of a considerable rate for the decay 
p -^ /i^7r° will be a clear indication for a gauge mediated proton decay. 

One can say in general, that the branching ratios of the nucleon decay can teach us about 
the "fundamental" mass matrices as they depend on all mixing angles. The present huge 
freedom in the mass matrices would then be strongly restricted and one could better 
understand the origin of the fermionic masses. 
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